arXiv:1509.06876v3 [astro-ph.HE] 4 Mar 2016 


Mon. Not. R. Astron. Soc. 000, ^{22] (2002) Printed 7 March 2016 (MN style file v2.2) 


Capturing the electromagnetic counterparts of binary 
neutron star mergers through low latency gravitational 
wave triggers 

Q. Chu* * * § ^, E. J. Howellf^, A. Rowlinson H. Gao^, B. Zhang®, S. J. Tingay^’® 
M. Boer^ and L. Wen^§ 

^School of Physics, University of Western Australia, Crawley WA 6009, Australia 
^ CSIRO Astronomy and Space Science, Sydney, Australia 

^ Anton Pannekoek Institute, University of Amsterdam, Postbus 94249, 1090 GE, Amsterdam 

'^The Netherlands Netherlands Institute for Radio Astronomy (ASTRON), PO Box 2, 7990 AA Dwingeloo, The Netherlands 
^Department of Astronomy, Beijing Normal University, Beijing 100875, China 
^Department of Physics and Astronomy, University of Nevada Las Vegas, NV 89154, USA 
"^International Centre for Radio Astronomy Research, Curtin University, Perth, WA 6845, Australia 
^ARC Centre of Excellence for All-sky Astrophysics (CAASTRO) 

^CNRS - ARTEMIS, boulevard de I’Observatoire, CS 34229, 06304 Nice Cedex 04, France 


7 March 2016 


ABSTRACT 

We investigate the prospects for joint low-latency gravitational wave (GW) detection 
and prompt electromagnetic (EM) follow-up observations of coalescing binary neutron 
stars (BNSs). For BNS mergers associated with short duration gamma-ray bursts 
(SGRBs), we for the first time evaluate the feasibility of rapid EM follow-ups to capture 
the prompt emission, early engine activity or reveal any potential by-products such 
as magnetars or fast radio bursts. To achieve our goal, we first simulate a population 
of coalescing BNSs using realistic distributions of source parameters and estimate the 
detectability and localisation efficiency at different times before merger. We then use 
a selection of facilities with GW follow-up agreements in place, from low-frequency 
radio to high energy y-ray to assess the prospects of prompt follow-up. We quantify our 
assessment using observational SGRB flux data extrapolated to be within the horizon 
distances of the advanced GW interferometric detectors LIGO and Virgo and to the 
prompt phase immediately following the binary merger. Our results illustrate that 
while challenging, breakthrough multi-messenger science is possible with EM follow¬ 
up facilities with fast responses and wide fields-of-view. We demonstrate that the 
opportunity to catch the prompt stage (< 5s) of SGRBs, can be enhanced by speeding 
up the detection pipelines of both GW observatories and EM follow-up facilities. We 
further show that the addition of an Australian instrument to the optimal detector 
network could possibly improve the angular resolution by a factor of two and thereby 
contribute significantly to GW-EM multi-messenger astronomy. 

Key words: binaries: close - gravitational waves - gamma-ray burst: general - 
methods: observational - stars: neutron 


1 INTRODUCTION 

The LIGO detection of GW150914 from the inspiral and 
merger of a p air of black holes has ignited a new era 
of astronomy ll Abbott et al.l [2016al l. To fully exploit this 
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new frontier it has been expected that triggered follow-ups 
by facilities operating outside the GW spectrum will be¬ 
come common place. During this new multi-messenger era, 
prompt electromagnetic follow-ups from low-latency GW 
triggers can provide a pan-spectral pathway of understand¬ 
ing into the inner engines and underlying mechanisms be¬ 
hind a multitude of GW sources. _ 

The Advanced LIGO (aLIGO; lAasi et ^ l2015ll detector 
launched its first science operations at sensitivities that 
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would allow the detection of coalescing systems of binary 
neutron stars (BNSs) in the range of 40 -80 Mpfl aLIG O will 
be followed by Advanced Virgo ('AdV: lAcernese et al.ll^isl 'l 
later this year (2016). By 2019 aLIGO is expected to reach 
design sensitivity, with a BNS detection range of 200Mpc. 
AdV is expected to reach this milestone by 2018 with BNS 
obser vations out to a distance of ISOMpc (lAcernese et ahl 
l2015h . Other second-generation detectors are also planned or 
undergoing construction. KAGRA, a Japa nese instrument , 
is envisioned to begin operation in 2017 iIAso et al.ll2013l) 
and fulfill its base line design arou nd 2018-2019 ( Somival 
l2012f) . LIGO-India jAasi et al.ll2013l ) has obtained approval 
and is expected to be operational from 2020, reaching a de¬ 
sign s ensitivity at the same level as aLIGO no earlier than 
2022 (lAasi et al.ll2013l) . 


Goalescing systems of BNSs are promising sources of 
GWs due to the fact that the inspiral phase is well mod¬ 
elled by the post-Newtonian (PN) approximation. The rel¬ 
atively low masses of BNS systems mean that these events 
will span most of the frequency band of aLIGO detectors. 
At design sensitivity, aLIGO is predicted to make 0.4 to 400 
BNS detections per year according to curren t rate predic¬ 
tions jAbadie et al.ll20ld : ICoward et al.lI 2 OI 2 I I . The coales¬ 
cence rate of of binary systems comprising of neutron stars 
and black holes (NS/BH) and binary black holes (BBH) are 
less confidently predicted due to a sparsity of observational 
data. Studies have predicted yearly detection rates of order 
0.2 - 300 detections for N S/BH events and 0.4 - 1000 for 
BBH events llAbadie et akll^lOl L For systems including at 
least one BH the detection of GW150914 is important as it 
proved that BH birth kicks do not always brea k up the evolu¬ 
tionar y formation channels for these sources (lAbbott et al.l 

l2016blL 


Multi-messenger astronomy is one of the highly prior¬ 
itized research areas in the GW community (I Howell et al.l 
l2015h . Gomplementary observations from electromagnetic 
(EM), neutrino or high energy particle facilities would not 
only improve the confidence of a GW detection, but also 
maximise the science return by providing a wealth of ad¬ 
ditional astrophysical information. A number of plausible 
EM emission mechanisms and end products that could be 
produced during the inspiral and merger phase of coalescing 
compact objects have been proposed. We summarise some of 
the popular scenarios in Figure [T] predictions are highly un¬ 
certain as illustrated by the various pathways and end prod¬ 
ucts and could involve combinations of BHs and/or magne- 
tar formation or Fast Radio Bursts (FRBs). The anticipated 
onset time, duration, and observable wavelengths shown sug¬ 
gest that transient EM emission could po ssibly occur within 
a few seconds after the b inary merger (iMetzger fc Berg^ 
I2OI2I : [Centrella et al.ll20l3 j. To catch these events, generat¬ 
ing GW event triggers without delay (thus low-latency) and 
obtaining prompt EM follow-ups are crucial. Achieving joint 
GW-EM observations of these events could reveal the major 
processes and interactions at play during the binary merger 
process (see detailed discussion in section O. 


^ This range refers to the detectable range averaged over all sky 
locations and source orientations for BNS systems with individ¬ 
ual masses of 1.4 Mq. It is approximately 1/2.26 of the horizon 
distance - the range for optimally oriented and located sources. 


A significant effort from both GW and EM astronomy 
community has been undertaken to a chieve prompt GW- 
triggered EM followup observations (lAbadie et al.l 1201 ll : 
IShawhanll2012ll . though as yet, no joint GW-EM detections 
have been made. Recent developments in the areas of astro¬ 
nomical and computing technologies in the following three 
fronts will further improve the chance of catching EM tran¬ 
sients during the binary coalescence: 

(i) Advances in low-latency GW search pipelines that are 
capable of generating event triggers within minutes upon the 
arrival of a detectable signal. 

(ii) Localisations are likely to improve as more GW de¬ 
tectors are constructed. 

(iii) Development of wide-held EM telescopes capable of 
rapid response within tens of seconds and with helds of view 
(FoV) comparable to the error area of GW source sky direc¬ 
tion (typically tens of square degrees for the aLIGO/AdV 
network). 


The hrst point has been studied since the initial LIGO 
era and a number of low-latency GW trig ger-generation 


pipel i nes have been proposed and tested (lAbadie et al 
20]J ; iBuskulic et al.ll2010l : I Gannon et al.l l20ld : iLuan et al 


2011 ; iHooper et al. 201lll : th e prospects of ra pid alerts from 


GW events are discussed in IShawhanI (|2012| L 

The benehts of a larger network given in the sec¬ 
ond point have prompted studies of source localization 
by advanc ed networks by a number of authors (see the 


2 OI 2 I: IVitale & ZanolirJl201ll: Schutd 201ll: iGhu et al. l2012|j. 

Wen et al.l ll2008f): Fairhurst (12009 

1: Wen & Chen 1 I 2 OIOII: 

Zanolin et al.l ll2009f): lAbadie et al. 

ll2010f): iNissanke et al.l 


of source localisation and estimation of localisation accuracy. 
For low-latency detections, studies of the early localization 
of BNS sources when their GWs are detected 10s of s econd s 
before the merger ca n be found in iManzotti fc Dietd ll2012ll : 
ICannon et al.l ll2012t) . 

_ A recent study of the sky localization (ISinger et al.l 

I 2 OI 4 II for the early detector network expected in 2016-17 in¬ 
dicate that the error in GW localization will be of order 200 
deg^. Methods of identifying GW sources from a large sky 
error regions will therefore have g r eat importance a n d hav e 
been proposed by Nissmtej^alJ ll2013ll ; iFan et aP ll2014h : 


ICowperthwaite fc Berger! ( 20151 1. 


To combat the large error regions the range of wide-field 
EM instruments that are in operation or are being developed 
across the entire range of the EM spectrum will be impor¬ 
tant. Additionally, the possibility of capturing the early EM 
emissions can be provided through low-latency alerts; this 
will require fast response EM instruments. A range of such 
instruments from high to low energy will be fully discussed 
in section 4. 

This paper will investigate the types of breakthrough 
science that will be possible through the low-latency de¬ 
tection and localisation of BNS GW events, in particular, 
how these will enable prompt GW-EM observations of BNS 
coalescence in the era of advanced detectors. While the sci¬ 
ence of GW-EM astronomy has been widely investigated 
previously for longer time sca le EM emissions ll Howell et al.l 
I2OI5I : IMetzger &: Bergeill2012l . e.g.) the feasibility to capture 
prompt and early EM emissions from binary mergers has not 
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Figure 1. A cartoon illustrating some of the possible multi-messenger pathways and end products of coalescing systems of NSs and 
BHs in different energy bands. We also show the approxima te timescales from GW merger. Short duration gamma ray bursts (SGRBs) 
are strongly linked with the merger of of compact objects l|Paczvnskiiri986l : lEichler et al.|[l989al : iNaravan et al.|[l992al : iRezzolla et al.l 
2Qlll:lGehrels et al.ll2Q05l: iBereer et al.|[20Q5l : iBloom et al.ll200i^ and could be accompanied by an optical flash or a fast radio burst (FRB 


Totanill2013h . If a stable or supra-massive magnetar is formed, the long lived X-ray plateaus ob served in many SGRBs could be the 
signature of direct dissipation of magnetar wind energy dRowfinsoi^^^ 2^^, 120151 : [Zhang|[2Ql3l : the collapse of a merger product to 


a black hole could also produce an FRB llZha nd 12014: iFalck e^^RezzoUa]^^^ . At later times predictions for an optical/IR kilonova 
llLi &: Paczvhsklll998l: iMetzger et al.ir2010l: lBarnes^rKasenll2Q13ll have gained support though recent observations llBerger et al.|[^13l: 
[TanvI^^aLll^l3[). For a relativel y small opacity, a kilonova precursor may appear early enough which would benefit from the low latency 
observations llMetzger et al.ll201fih . 


been studied before. To examine the possible EM responses, 
we use a comprehensive simulation that extends on previ¬ 
ous investigations of the early detection and localization of 
compact binary coalescence (CBC) sources by extensively 
sampling a range of parameters of the NS waveform and 
sky localisation. We produce statistics showing the expected 
number of detections and their sky localisation error regions 
at different times before merger. Based on these results, we 
for the first time extrapolate the temporally complete EM 
data currently available for SGRBs to investigate the ca¬ 
pacity of EM instruments for short time scale EM follow-up 
observations through the 7 -ray, X-ray, optical, and radio. We 
consider a range of possible emission mechanisms for BNS 
mergers; some of the possible pathways and end products are 
illustrated in Figure [T] Furthermore, our investigation uses 
the full compliment of GW detector network configurations 
that will be available through the advanced detector era by 
including KAGRA, LIGO-India and a proposed southern 
Hemisphere detector AIGO. 

The paper is organized as follows: in section [2] we de¬ 
scribe the framework we use to simulate the GW detection 
and localization and provide these results in section [S] In 
section[4]we discuss the potential of EM telescopes for rapid 
follow up observations and in section [S] provide a review 
of some of the possible EM emissions that could be cap¬ 
tured by prompt observations. Then in section [ 6 ] we use our 


simulated detection and localisation statistics with observa¬ 
tional data to determine the astrophysical implications of 
low-latency GW detection and prompt EM follow up. We 
then summarise our findings in section [7] 
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Table 1. The workflow of a low-latency CBC detection pipeline showing the latencies achieved in science and engineering runs as well 
as a projected optimistic latency. The LIGO S6-Virgo VSR2/VSR3 joint observing run was performed during 2009-2010; ER5 and ER6 
were engineering runs used to streamline data analysis pipelines by simulating the environment of a real science run. | This assumes that 
vetting can be eliminated through automatic processing. 


Process 

Latency in S6-VSR2/VSR3 joint run 

ER5/ER6 t latency 

Optimistic latency 

GW data acquisition, calibration Sz distribution 

1 min 

'^lOs 

0(l)s 

GW trigger generation 

2>^7 mins 

'^40s 

0{l)s 

Follow-up preparison 

2'^3 mins 

60s 

0(l)s 

Human vetting 

20~30 mins 

N/A 

Os t 


2 METHODOLOGY 

2.1 GW Detection 

The response of a laser inteferometric GW detector to a GW 
signal at a given time t is given by: 

h(t) = F+{t)h+{t) + (1) 

where F+ and Fx are the antenna pattern functions of two 
independent plus and cross GW polarizations denoted as h+ 
and hx- We adopt the conventio n of the antenna responses 
given in I.Taranowski et all (|l998l 'l . The antenna patterns are 
expressed as functions of the source sky location—the right 
ascension a, the declination 5 and the polarization angle ip. 
In our simulation, we calculate the waveform h+{t),hx{t) 
using the second order post newtonian (PN) expansion. It is 
a function of eleven parameters: the luminosity distance of 
the source r, the orbital inclination angle t, the coalescence 
phase (pc, masses and spins (three parameters for each spin) 
of these two inspiral bodies. In this paper we only consider 
non-spinning BNS coalescence. 

Matched templates are constructed from the CBC pa¬ 
rameter space to perform hltering o n the detec tor data to 
search for and extract GW signals (|FinrJfl992l 'l. The out¬ 
put of matched filtering is the signal-to-noise ratio (SNR), 
from which a selection criteria is applied to yield credible 
GW triggers. Throughout this paper we assume that our 
detections yield an optimal SNR for each source; that is we 
assume the template matches the source exactly. 

For a network of detectors, two strategies can be con¬ 
sidered to select GW candidate triggers. The coincident se¬ 
lection criteria (which we refer to as the triggering criteria 
in the later part of the paper) requires at least two detec¬ 
tors to detect signals above a predefined SNR threshold. A 
coherent selection criteria assumes that data from different 
detectors form coherent responses from the source; thus a 
single statistic is produced and applied against a threshold 
for trigger generation. 

2.2 Latency of GW Event Triggers 

The latency of a GW detection pipeline is defined as the time 
elapsed from when the source has sufficient SNR for detec¬ 
tion to the time it is registered in the GW event databas^l. 
Table [T] outlines the workflow of a typical low-latency GBC 
pipeline showing the latencies achieved in the LIGO S6- 
Virgo VSR2/VSR3 joint observing run, the ER5 and ERG 
engineering runs as well as a projected optimistic pipeline. 

^ GraceDB; The Gravitational-wave Candidate Event Database 


The workflow was fi rst practiced in th e latest joint 
LIGO/Virgo science run jAbadie et al.l I2OIIII and the la¬ 
tency is shown in the second column of Table [T] In this 
process detector data was acquired, conditioned, and dis¬ 
tributed. Gandidate triggers including CBG triggers were 
generated by different low latency GW search pipelines and 
submitted to the database. Significant events were then se¬ 
lected and validated manually to form alerts to participating 
EM observatories. The total time lag was of order 30 mins of 
which 20 mins consisted of human validation. The third col¬ 
umn of Table [1] shows the results of a recent engineering run 
with low-latency pipelines specific to search for GBC signals 
and the final column shows estimated optimistic latencies 
possible in the future. 

GW detector data are first acquired onsite and cali¬ 
brated and distributed to the computing node where search 
pipelines are running. From the latest engineering runs, this 
process takes of order 10 seconds. In principle, there is no 
intrinsic obstacle preventing this procedure to achieve min¬ 
imal latency. Therefore, we assume in the most optimistic 
scenario that data can be packed, calibrated and distributed 
within 0(1) s. 


At the time of writing, three low-latency GW search 
pipelines have been implemented for t he LIGO-Virgo col- 


laboration to search for GBC signals llAbadie et al. 

2012; 

Buskulic et al.l 2010l: ICannon et al.ll2012l: Luan et al 

2011; 

Hooner et al.l l201lfl. These three search Dinelines 

have 


all achieved an average of around 40s latency from the 
time of data acquisition to trigger registration. This la¬ 
tency is mainly due to the processes of data whiten¬ 
ing, filtering, and coincident analysis. By employing high- 
performance c omputing hardware such as graphic p rocessing 
units (GPUs) (iGhung et al.|[201ol : iLiu et al.ll2012l l and new 
high-performance algorithms, this latency could be further 
reduced. We assume in the future 0(l)s latency for filter¬ 
ing and trigger generation from multiple detectors can be 
achieved in the most optimistic case. 

The third row of Table [T] shows the latencies introduced 
through follow-up preparation, which includes the gather¬ 
ing of event parameters, e.g. sky localization and calculat¬ 
ing telescope tilings. The main latency of this part comes 
from event localisation. A low latency local ization pipeline 
BAYS TAI0 has been developed and tested (ISinger fc Pried 
I2OI5II and has achieved latencies in localization of seconds; 
the run time of this method is shown to be almost perfectly 
proportional to the number of GPU cores. Thus this latency 
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Figure 2. The parameter distributions for the detected sample of sources. We show the distributions of: (a) sky location;(b) inclination; (c) 
chirp mass and (d) total mass. 


could be largely reduced and in the optimistic scenario we 
assume 0(l)s latency for this procedure. 

Human vetting introduces the highest latency as the 
first detections needs to be validated carefully. It is nat¬ 
ural that once enough detections have been recorded, the 
knowledge gained will allow this part of the routine to be 
eliminated by vetoing candidates automatically. We there¬ 
fore estimate 0(1) s for the optimistic case. 

In summary, we consider several latency scenarios that 
include both GW detection and localization for GW trigger¬ 
ing: an ideal zero latency and latencies of [Is,10s,40s] which 
takes into account the range of extreme to plausible laten¬ 
cies to send out a GW trigger. In section [5l our discussion is 
focused on a latency of 40s which is achievable in the coming 
science runs and an optimistic latency of Is. 


2.3 Sky localisation accuracy 

To determine the localisation acc uracy of the detected GW 
events, we utilize the formalism of lWen fc GhenI (l20inl l who 
used the Fisher information matrix to derive a geometrical 
expression for sky l ocalization acc uracy. According to the 
Cramer-Rao bound (|Cramerlll946h . the lower bound of the 
variance of any unbiased estimator is determined by the in¬ 
verse of the Fisher matrix F. Using this property, one can 
calculate the lower bound on the localisation error area. This 
estimation is applicable to any unbiased estimators and to 
all GW waveforms. This method is also used in ISchut j 


ll201lll for a sky localization comparison of different detec¬ 
tor networks. _ 

IWen fc ChenI ll201(]ll presented angular resolution for¬ 
mulae for two different detection scenarios: a best case sce¬ 
nario assumed that the waveform was known but the ar¬ 
rival time was unknown; a worst case scenario assumed that 
both these quantities were unknown. The calculations in¬ 
clude both phase and amplitude information of the signal. 
The angular resolutions derived from the best case scenarios 
were shown to be around half the areas of those de termined 
assuming the worst case scenario jChu et al.ll2012h . In this 
paper, we assume the best case scenario to calculate angular 
resolution for each detection. This is given as: 


AH = 4v/27i 


(E.6 


E 


J,K,L,M 


X vml) ■ 


( 2 ) 


with n the direction of the source which is a function 
of right ascension a and declination 5 and rxj the dis¬ 
tance vector from detector K to detector J; the quantity 
= /Eo interpreted as noise-weighted en¬ 

ergy flux received by the ith detector and z is the SNR (as 
discussed in section EHJ- 


2.4 Monte-Carlo Simulation 

To produce a sample of detection statistics we simulate a 
population of 20,000 BNS sources. Our Monte-carlo code 
extends other studies by considering a full complement of 
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parameters for the BNS waveform and the sky location. The 
values of right ascension a and polarization angle ip are draw 
from a uniform distribution within the range [0, 2-k]. Values 
for the declination <5 angle and inclination angle l are drawn 
from a distribution of cosine values uniformly distributed 
over the range [—1,1], For each event a source distance r 
is chosen from within [10, 500] Mpc; this range encompasses 
the horizon distances of 428Mpc for aLIGO, 332Mpc for 
AdV, and 322Mpc for KAGRA. For the component masses 
of each event we assume a Gau ssian distribution based on 
the study of lOzel et ^ (120121 '): this study suggests that 
the mass distribution of BNS systems peaks around 1.33M0 
with a standard deviation of O.OSMq. 

It is expected that the aLIGO-Livingston (denoted L), 
aLIGO-Hanford (H), and AdV (V) detectors will reach 
their full sensitivities around 2019. To model the sensitiv¬ 
ities of the advanced LHV detector network, we use the 
‘zero tuning, high power’ spectral density noise curve of 
ILIGO Scientific CollaborationI (l2009l 'l for aLIGO-Livingston 
and aLI GO-Hanford; for AdV we use the spectral density 
curve of IPunturol ll2012f) . 

Our triggering selection criteria requires a SNR of 7.5 
in at least 2 detectors 0- This is equivalent to 91% of the 
reach of a single aLIGO detector triggering off an SNR of 
8; this also translates to a 194 Mpc horizon distance. As¬ 
suming there will be tens of detections of an LHV network 
per year, we present our detection estimates as percentages - 
this allows the statistics for different thresholds to be easily 
obtained through scaling. 

We obtained a total of 1010 LHV coincident detections 
from our simulated population of 20,000 sources. For an 
enlarged detector network we assume the same coincident 
threshold. The key statistics will be presented in section 
We note that some studies have used a coher ent analysis 
threshold to select sources, e.g. ISchut j |[20lJ). The coin¬ 
cident selection thresholds shown here can be simply con¬ 
verted to corresponding coherent analysis thresholds via the 
detection volume. The statistics of detections from coher¬ 
ent analysis should be similar to the statistics of detections 
shown in the following section. 


3 RESULTS FOR GW TRIGGER EVENTS 

3.1 Parameter Distributions of Detected GW 
Events from LHV network 


Fig.2 shows the distributions for the different parameters 
of the 1010 simulated BNS detections for an LHV network 
in the advanced era. Fig l2al shows that the distribution of 
detected sky locations is bimodal. The locations of max¬ 
ima and minima are consistent with the landscape of the 
antenna power pattern of a LHV detector network (e.g. 
Fig. 2 in Nissanke et all 2009ll. As shown i n several stud- 
ies I Manzotti fc Dietd I2OI2I : ISchutd I2OIII : iNissanke et al.l 
12001 ) adding more detectors to GW search will significantly 
improve sky coverage. 

The distribution of detected inclinations shown in 
Fig l2bl is bimodal and shows a detection preference towards 


We note here that the detection pipeline gstlal is configured 
to use a single detector threshold SNR of 4. 


low inclination angles. We note however that the fraction of 
dete cted inclinati ons below 25° is less than one quarter (see 
also ISchut jl201ll L The distributions of the chirp mass and 
total mass of sources are shown in Fig. [2^ and Fig. I2dl 


3.2 The localisation of BNS detections 


In this section we present the localization evolution statis¬ 
tics for the LHV detections and for an extended detector 
network that covers the entire parameter space of possible 
BNS sources. There is a crucial interplay between sending 
out a prompt alert to an EM facility and surveying the GW 
error region before the EM source fades. An early trigger can 
allow a telescope to get on source, but at the expense of a 
larger GW error region. The situation will vary significantly 
for different types of telescope and is also dependent of the 
type of EM emission. In preparation for the later analysis of 
EM follow-up in Section 6 we have calculated the percentage 
of sources detectible within different sized error regions at 
different times before merger. We note that previous studies 
have calculated the evolution of the localiz ation property 
with time using a qualitative analysis (e.g. ICannon et al.l 
I2OI2II rather than the quantitative approach adopted in this 
paper. 

Figure |3] displays the percentile charts for the evolution 
of network SNR (I3all and localisation error area (I3bl) for the 
1010 LHV detections as they approach merger (t=0). The 
network SNR is calculated as root-sum-square of individual 
SNRs. The plot shows that the median network SNR (thick 
solid line) is less than 5 at 100s before coalescence. This SNR 
increases by 1 every 10 seconds during -100s to -10s and by 
0.5 every second during -10s to -Is before reaching a SNR of 
16 at merger. We find that 10% of the detections will yield 
a SNR greater than 27 at merger time which corresponds to 
a localisation area of around 5 deg^. 

Figure [3b] shows the percentile range of localisation er¬ 
ror areas at 90% confidence level as a function of time be¬ 
fore merger. On average, the angular resolution is about 
30000 deg^ at -100s. As the network SNR increases, this im¬ 
proves to < 5000 deg^ at -40s, < 500 deg^ before -10s and 
less than a hundred square degrees at -Is. These results sup¬ 
port the early locali sation results of a number of recent stud¬ 
ies: se e for examp le lCannon ef^aJjjr^Oljl. [Manzotti fc Dietj 
1 I 2 OI 2 II and Ch.3, Fig 3.3 of lSin^ (l2015h . 


3.3 Improved Early Warning Events with a 
Larger Detector Network 

In this section we investigate the benefits of an enlarged 
GW network in terms of early detection and localisation; 
the scenarios considered are applicable to the advanced de¬ 
tector era beyond 2019. We consider a number of additions 
to the LHV network: these are KAGRA (labeled as J), the 
proposed Indian detector: LIGO-India (labeled as I), and 
AIGO (labeled as A). Though AIGO is not approved finan¬ 
cially at the moment, it is considered here for comparison. As 
in section [ 2 ^ we use the aLIGO spectral density for LIGO- 
India and AIGO, and the spectral density from lAso et ^ 
1 I 2 OI 3 II for KAGRA. The main results are given in Table [5] 
and Figure [4] 

The second column of Table [2] shows the detection ca- 
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Table 2. The detection statistics obtained for different advanced detector networks shown in the first column. The second column shows 
the relative detection capability of each network in comparison to LHV. Columns 4 to 7 provide two sets of statistics at different times 
before merger for each detector network: in the first row we provide the percentage of detections achieved pre-merger; the second row 
gives the sizes of the error regions required to contain 10%, 50% and 90% of the detections at exact times before merger. | This number 
is lower than that of LHVI/LHVA network as the sensitivity we adopted for detector J is 0.75 of that of detector I/A. 


Network 

Detection 

capability 


40s 

Time before 

10s 

merger 

Is 

Os 

LHV 

1.00 

Percentage detected 

9% 

35% 

76% 

100% 



[10%/50%/90%] error regions [deg^] 

300/1000/2200 

70/270/780 

20/80/240 

5/20/50 

LHVJ 

1.12 t 

Percentage detected 

8% 

32% 

76% 

100% 



[10%/50%/90%] error regions [deg^] 

240/820/1600 

50/200/470 

10/40/120 

3/10/30 

LHVI 

1.28 

Percentage detected 

9% 

33% 

75% 

100% 



[10%/50%/90%] error regions [deg^] 

120/380/910 

40/130/350 

10/40/100 

3/9/20 

LHVA 

1.27 

Percentage detected 

9% 

35% 

76% 

100% 



[10%/50%/90%] error regions [deg^] 

80/310/850 

30/90/250 

9/30/80 

2/6/20 

LHVJI 

1.35 

Percentage detected 

8% 

32% 

75% 

100% 



[10%/50%/90%] error regions [deg^] 

100/340/830 

30/110/300 

9/30/80 

2/7/20 

LHVJIA 

1.57 

Percentage detected 

9% 

33% 

75% 

100% 



[10%/50%/90%] error regions [deg^] 

60/180/400 

20/60/140 

6/20/40 

1/4/10 


pability of different detector networks in terms of the num¬ 
ber of detections normalized to that obtained from a LHV 
network (the reference detection capability of LHV is given 
here as 1.00). For larger network configurations, we find that 
the detection capability is improved considerably. Using the 
same coincident triggering threshold for different detector 
networks, we obtain 1133 detections for LHVJ, 1291 detec¬ 
tions for LHVI, 1366 detections for LHVJI, 1280 detections 
for LHVA, and 1584 detections for LHVJIA. The detection 
capability of a larger detector network depends on the sen¬ 
sitivities of the additional detectors. We find that adding a 
single detector with the same detection capability as aLIGO 
to the LHV network will increase the detection rate by about 
30% while adding KAGRA to the LHV network will improve 
the detection rate by 12% as the sensitivity of KAGRA is 
0.75 of that of aLIGO. The detection percentage for each 
network at different times of the final merger phase is shown 
in the hrst row for each network configuration. 

Table[2]also shows localisation achievable at exact times 
before merger for different GW network configurations. We 
illustrate this by showing the sizes of the error regions re¬ 
quired to contain 10%, 50% and 90% of the detections. All 
the error region numbers are estimated at the 90% confi¬ 
dential level. We find that adding a KAGRA or LIGO-India 
instrument to the LHV network can improve the localisation 
by a factor of 2. Increasing the network to five detectors can 
lead to a factor of 3 improvement. Adding AIGO on its own 
results in an improvement that is equivalent to adding both 
KAGRA and LIGO-India due to the longer baseline and op¬ 
timal location. Adding AIGO to the 5-detector network, can 
further improve the localisation by a factor of around 2 - this 
is a factor 5 improvement over the initial LHV network. 

An important fact is that only a fraction of detections 
can be detected at times prior to merger. We take this into 
consideration and show the localisation error area for the de¬ 
tections at different times prior to merger. Fluctuations of lo¬ 


calisation error area at early times in Fig. U are due to small 
number of detections in that epoch. We note our estimated 
localisation error areas determined near merger time are 
compa rable, though sli g htly m ore optimistic than iFairhursd 
(l2010ll and lAasi et ^ (l2013l L Whilst the aforementioned 
studies used the Fisher matrix for the triangulation localisa¬ 
tion method, our estimates are based on the Fisher matrix 
regardless of the localisation method, thus represent a lower 
limit of all non-biased estimations. 

In practical terms, the latencies introduced in the GW- 
EM effort will affect the timing of the EM followup obser¬ 
vation and therefore constrain the phenomena we are able 
to observe. Using Tabled one can easily consider different 
GW-EM reaction latency scenarios and find the detection 
ratio and localisation error area for follow-up observations 
at merger or at post-merger timescales. For example, as¬ 
suming a 40 second total reaction latency - incorporating 
the whole GW trigger generation and EM-followup readiness 
latencies - can be achieved, those GW events that pass the 
GW triggering criteria 40s before the merger can potentially 
be observed right at their merger time by EM instruments 
with wide FoVs and rapid responses. For the LHV network, 
which is expected to be operational at early stage of ad¬ 
vanced era, around 9% of all detections can be observed at 
the merger time within 40s. The median localization error 
area of the detections at the merger time is over 1000 deg^. 
The median localisation error area can be brought down 
to about 400 deg^ by an LHVI detector network and 10% 
detections of LHVI network at the merger time could be lo¬ 
calised within 120 deg^. With a five-detector network, 10% 
of the detections could be discovered and localised within 
100 deg^ at the time of merger. 

If a greater effort could reduce the total reaction la¬ 
tency to 10 seconds (this highly optimistic scenario would 
require an EM facility to be on source), the detection ratio 
at the merger time will be raised to around 33% for the six 
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(a) 



(b) 

Figure 3. Percentile curves of (a) network SNR and (b) Af2 (er¬ 
ror box at 90% confidential level) as a function of time before 
merger for the simulated BNS triggers obtained by an LHV net¬ 
work. 



Figure 4. Median localisation error areas of detections at exact 
times prior to merger for different detector networks. The fluc¬ 
tuations seen at early times are due to the smaller number of 
detections in that regime. 


networks which is three times the observed detections at 40 
second latency. More importantly, the localisation accuracy 
is also three times better than that achieved with 40 sec¬ 
ond latency. For the LHV network, the median localization 
area is improved 4-fold to 270 square degrees and 10% of 
detections could be localised within 70 deg^. Adding KA- 
GRA or LIGO-India into the LHV network will improve the 
localisation to less than 50 deg^ for 10% of their detections. 
The localisation accuracy will be further brought down to 
about 30 deg^ for a network of LHVA or LHVJI of 10% of 
their detections. The best localisation accuracy comes at an 
operation of six detectors with around 90% of the detec¬ 
tions localised to within 140 deg^. We consider further the 
possibilities of existing EM facilities in the next section. 


4 EM TELESCOPES FOR PROMPT GW 
FOLLOW-UP OBSERVATIONS 

In this section we will discuss some of the EM instruments 
in different EM wavelength bands that will be able to facil¬ 
itate prompt GW follow-up - we consider a selection from 
those that currently have MoUs for aLIGO/AdV follow-ups 
in place. We focus first on the high energy which is tradition¬ 
ally low latency (usually employing triggering instruments) 
but will not be fast enough unless there is an agreement for 
prompt Target of Opportunity (ToO) in place; we note that 
high-energy observations of GRBs over the last two decades 
have produced significant science returns in this domain. We 
then discuss optical telescopes which is proven for low la¬ 
tency through the successful use of fast response instruments 
for GRB follow-up. Finally we consider radio which though 
electronic steering and/or rapid slew-times, can make rapid 
response follow-ups achievable. 


4.1 X-ray and 7 -ray follow-ups 

Table [3] shows the various properties of some of the high 
energy instruments with LIGO/AdV MoUs in place. 

Space based 7 -ray follow-ups will be led by Fermi and 
Swift. These instruments provide the capability for the dis¬ 
covery of simultaneous EM counterparts; in this respect the 
large FoV of the Fermi-GBM, which sees almost the entire 
unocculted sky can be advantageous. A typical follow-up 
procedure will be to search through the data for a coincident 
signal from a location of the sky consistent with the LIGO 
error region. Assuming coincident [GW/Fermi onboard] trig¬ 
gers, the response time will be around a few minutes. If 
there are no on-board triggers, a sensitive search through 
the Fermi data set will require of order 4-12 hours. The 
Fermi-LAT can provide more information of the high en¬ 
ergy afterglow. Typically the LAT autonomously re-points 
after a GBM trigger with a processing delay of around 30s. 
If no on-board trigger is available, a ground based search 
of the data can be initiated; the latency will depend on the 
time between the GW candidate and the data downlinl0 
and processing, of order one to a few hours. 

The Swift 7 -ray instrument is a wide-field (1.4sr) coded 


® The communication going from the satellite to the ground 
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Table 3. The prope rties of some of the hi^h -ene r^y instrum e nts pl anned to be operational during aLIGO. fll iGehrels et all l|2Q04^ : [2] 
iMeegan et alj l2009lb [3l iBartos et alJ ll201^') : [4l iGdtz et alJ ||201^') . NOTES: § Sensitivity approximated using IZh an d J2013t) equation 

1 3) -Fy th = 6 X 10“^^/Tobs erg sec~^ cm~^and assuming T^hs = 10s. | Sensitivity over 100s estimated using Fig 2. of I Atwood et all 
2003) and conve r ted to ergsec”^ cm~^ using a standard band function with a = —1, ^ = —2.25 and Ep = 511 keV. b Approximated 
using [plmk et ahl ll2013l) . t] Sensitivity in survey mode based on [3] in which exposure time includes an estimate of the required slewing 
times to tile a 1000 deg^ area using convergent pointing mode. 


Instrument 

Energy Range 

FoV 

[deg2] 

Sensitivity 
[ergsec”^ cm~^] 

Exposure Time 

Response Time 

Ref 

Swift - BAT 

15-300 keV 

4600 

1.2 X 10 "® 

_ 

Coincident Observation 

[1] 

Swift - XRT § 

0.3-10 keV 

0.15 

6 X 10-12 

10 

1-2 hours 

[1] 

Fermi - GBM 

8 keV -40 MeV 

30000 

4 X 10-® 

- 

Coincident Observation 

[2] 

Fermi - LATf 

0.02-300 GeV 

8000 

1.4 X lO-l" 

lOOs 

30s post GBM trigger 

[2] 

GTA b 

0.03 - 100 TeV 

6-8 deg^ 

6 X 10“^® 25GeV 

1000s 

20-60secs 

[3] 

CTA (Survey mode) 

0.03 - 100 TeV 

~1000 deg^ 

6 X 10-®@ 25GeV 

1000s 

20-60secs 

[3] 

H.E.S.S. b 

0.05 - 20 TeV 

15 

~ 6 X 10“® @ 25GeV 

1000s 

30s 

[3] 

SVOM - EGLAIRS 

4-250 keV 

89 

7.2 X 10-1° 

1000s 

~ 1-2 hours 

[4] 

SVOM - MXT 

0.2-10 keV 

64 arcmin 

5.6 X 10-11 

10s 

~ 1-2 hours 

[4] 


Table 4. The properties of a selection of the optical telescopes planned to be operational durin; 


Keller et al, 


Smith et al 


2007): [3l iGoward et ahl ll2010h: l4l|Hod a >p et ^ ||2004) : [5l ISteeghs fc Gallowavl ll201 
2014ldRellml ||2014|~) fSl lPaul et alj 1 I 2 OI 1 I) 


aLIGO_[l]jKl£tz_et_ay 1^0^ ); 
[6] iGhosh fc^elemana |201£ ); 


[ 2 ] 

[7] 


Telescope 

FoV 

[deg2] 

Limiting Magnitude 
R-band 

Exposure 

Ref 

TAROT 

3.5 

18 

60 

[1] 

SkyMapper 

5.7 

21 

110s 

[2] 

Zadko 

0.15 

21 

180s 

[3] 

Pan-STARRS 

7.0 

24 

30s 

[4] 

GOTO 

18.0-36.0 

21 

5m 

[5] 

BlackGEM 

40.0 

22 

5m 

[6] 

Zwicky Transient Facility (ZTF) 

47 

20.5-21 

30s 

[7] 

Ground based Wide-Angle Camera (GWAC) 

8000 

16 

10s 

[8] 


mask hard X-ray imager (Burst-Alert-Telescope; BAT). Al¬ 
though the Swift-^AT is more sensitive than the Fermi- 
GBM, it observes softer 7-rays; this last point is important 
for Fermi-GBM, follow-ups as SGRBs associated with BNS 
mergers have harder spectra than LGRBs. 

In terms of identifying transients, follow-ups in the X- 
ray have the advantage of a less crowded sky than that of 
the optical. For X-ray follow-ups the Sm/t-XRTs rapid slew 
time is particular advantageous; this has resulted in this 
instrument accumulating around 200 GRB redshifts though 
co-ordinated observations with ground based telescopes. The 
S'rui/t-XRT is able to respond rapidly to a BAT detection of 
new transient, localising to within 1-3 arcminutes and then 
distributing the position via the Gamma-ray Coordinates 
Network (GCN) within a 20 second window. 

Swift demonstrated the feasibility of rapid high energy 
ToO follow-ups by conducting the first multi-wavelength 
follow-up observation durin g the 2009-2010 LIGO/Virgo sci¬ 
ence run d Evans et al ][2 m3). Swift followed-up and observed 
the sky locations for two candidate GW events within 12 hrs, 
using a specially designed strategy to target the most proba¬ 
ble positions that the GW event originated from. To improve 
ToO requests with large positional uncertainties. Swift im¬ 
plemented tiled XRT observations in late 2011. However for 
this instrument, at the present time, a response time to ToO 


of order 1-2 hours (iKanner et al. I I2OI2II by the on board X- 
ray telescope (XRT) effectively rules this instrument out in 
terms of low latency follow-ups. If however, in the future, 
aLIGO could send triggered alerts directly to Swift without 
human intervention, a slew of order 10s of seconds would re¬ 
sult in early follow-up observations of the X-ray light curve. 

In addition to the Swift XRT, Astrosal[j, an Indian 
multi-wavelength satellite is was successfully launched in 
September 2015. Astrosat has instruments covering the UV 
(UVIT), Optical, Soft & Hard X-ray (the Cadmium Zinc 
Telluride Imager; CZTI). For follow-up operations, once a 
counterpart is identified, re-pointing will bring the target 
within the FoV of the other instruments for simultaneous 
coverage from optical though to hard X-ray. For low-latency 
follow-ups, CZTI, the hard X-ray coded aperture mask in¬ 
strument with a 36 deg^ FoV will be the most useful in¬ 
strument. As contact with the spacecraft is designed to take 
place once in an orbital period of 100 minutes, this is typi¬ 
cally is the fastest time scale over which a slew can be com¬ 
manded; ToO decision making may also be added to this 
response time. However, one low-latency scenario exists, in 
that if there was significant hard X-ray/7-ray flux in the 


http://astrosat.iucaa.in/ 
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Tabl e 5. The capacity of ra dio telescopes in consideration [Ij iTingav et alJ l|2013l '): pj lMurphv et aO l|2013h : [3] Ivan Haarlem et alJ i|2013l ') 
f; [4l lOosterloo et alJ 1 I 2 OIOII : t] APERTIF will be able to cover 100 deg^ at lower sensitivity in a “fly’s eye” mode. 


Name 

Frequency Range 

FOV (sq-degs) 

[sq-degs] 

limiting flux density 

Response Time 

Ref 

MWA 

80-300 MHz 

610deg2@150MHz 

20 mjy 

< 10s of secs 

[1] 

ASKAR 

700 MHz-1.8 GHz 

30deg^@1.4GHz 

640/iJy beam^ (10s int) 

< mins 

[2] 

LOFAR LBA (Inner) 

10-90 MHz 

450 deg2@60MHz 

5 mJy @60 MHz 

< min 

[3] 

LOFAR HBA (Core) 

110-250 MHz 

48 deg2@180MHz 

0.6 mJy @180 MHz 

< min 

[3] 

APERTIF 

1-1.7 GHz 

8deg^@lGHz (>100deg^)t] 

0.1 ^Jy 


[4] 


LOFAR sensitivities quoted for an observation comprising of a total bandwidth of 7.8 MHz and a 1 hour integration time assuming 6 
sub-array pointings are created, calculated using http://www.astron.nl/~heald/t6st/sens2.php 


(100-1000) keV band, the CZTI would detect the flash with¬ 
out re-pointing. Over this energy range the CZTI acts like 
a open detector with a 27r FOV; we note however that there 
will be no localization capability at these energies. 

Ground based Gamma-ray detectors can also contribute 
to the GW fo llow-up effort throu gh the Cherenkov Telescope 
Array fCTA: lDubus et al.ll2013ll . The CTA is a next gener¬ 
ation ground-based instrument that will begin construction 
around 2016, with full operation expected by 2020. CTA 
will improve over previous experiments ( H.E.S.S, VERI- 
TA^ and MAGICjj) through increased sensitivity, angular 
resolution (~ 2 at TeV), wider energy coverage (~ SOGeV 
- SOOTeV) and a larger field of view (6°-8°). This project 
will consist of two arrays: a southern hemispheric array fo¬ 
cusing on Galactic s ources and a northe rn hemispheric ar¬ 
ray on extragalactic. iBartos et al.l (l2014h has recently high¬ 
lighted some of the properties that make this instrument 
ideal for GW follow-up. These include the capability to ob¬ 
serve in a pointed mode with a relatively small FoV and in 
survey mode covering a larger area large sky area (~ 1000 
deg^ ); fast repo sitioning is of the order of a few seconds 
llDoro et al.ll2013l l. For low-latency follow-ups the ability to 
respond to ToO requests and st art monitoring the GW error 
ellipse within 30s is significant (iDubus et al. l2013ll . 

CTA is designed to have three types of telescopes with 
different mirror sizes to cover the full energy range - large 
scale telescopes (24m) covering < 300 GeV, medium scale 
(10-12m) covering 100 GeV-TeV and small scale (4-6m) cov¬ 
ering higher energies >10 Tev. Typically, extragalactic back¬ 
ground light (EBL), which is all the accumulated radiation 
in the Universe due to star formation processes, will atten¬ 
uate very high energy photons coming from cosmological 
distances through pair production with EBL photons. This 
produces the effect in that only the large scale telescopes 
will be triggered. However, for the case of a GW source 
withi n 200 Mpc, TeV pho tons are not expected to be annihi¬ 
lated llBartos et al.ll2014l ) ; thus all 3 types of GTA telescopes 
GTA could be triggered (lO’Brien fc Hinto^l2015h providing 
a greater energy coverage. 

H.E.S.S0, an array of 5 Gherenkov telescopes (with 4x 
12m and one 28m diameter mirrors) located in Namibia, 
has been operational since 2004. This telescope has a mean 


^ http://veritas.sao.arizona.edu/ 

® http://magic.iiippmu.mpg.de/ 

® http://www.mpi-hd.mpg.de/hfm/HESS/ 


time to get on sour ce from a random ob servation position 
of about 30 seconds llLennarz et al JliH). The capability of 
H.E.S.S. for GW follow-ups has been proven through prompt 
observations of GRB triggers. A follow-u p observation of 
GRB 070621 was achieved within 7 mins IIAharonian et all 
I 2 OO 9 II and GRB 100621A was observed within 10 mins. 

4.2 Optical follow-ups 

Optical follow ups have been hugely successfully in the last 
decade to discover the afterglows and host galaxies of GRBs. 
For Swift follow-ups, this strategy has relied on the prompt 
localisations of the XRT and the rapid distribution of coor¬ 
dinates via the GGN network. Recently, robotic optical tele¬ 
scopes have been able to identify the optical afterglows of 
Fermi GBM triggers that have large positiona l uncertainties 
(e.g. GRB 130702A: ISinger et aI.ll 2 oT I2OI5I I. A successful 
low-latency EM follow up to a GRB could potentially cap¬ 
ture an optical flash or the earliest stages of a particularly 
bright rising afterglow. 

Table 0] lists the capabilities of a number of optical in¬ 
struments that will be available during the era of advanced 
GW detectors. Gapturing an EM counterpart with such a 
fast response will be dependent on a number of factors re¬ 
lated to the telescope. The FoV of the instrument will deter¬ 
mine how quickly the GW localisation area can be surveyed 
(i.e. the number of tiles required); any subsequent detection 
will then be dependent on the limiting magnitude or sensi¬ 
tivity of the instrument and on the required exposure time 
to obtain the required s ensitivity (for an in depth discussion 
see iGoward et al.l[20l3 '). 

For prompt follow-up, the slew time of the instruments 
will be highly important; for telescopes with lower sensitiv¬ 
ities a rapid slew time can compensate through the ability 
to achieve pointing during the brightest part of the event. 
However, we note here that for optical instruments, a fast 
slew time is usually associated with a small FoV. 

Fast response optical telescopes such as TAROT have 
the ability to slew to a target position within 10 seconds. 
However, the FoV ~ 3.5 deg^ and exposure time of around 
60 s mean that time must be allowed to image the full error 
region. Assuming a 60 deg^ error box, corresponding to that 
of a LHVJIA net work at 50% confiden ce, using the frame¬ 
work presented in iGoward et al.l (l2014l L we find the imag¬ 
ing time will be of order 17 minutes. Larger FoV telescopes 
would be able to minimise the required imaging time. For 
example, Pan-STARRS, a system of four automated 1.8m 
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optical telescopes based in Hawaii has a FoV of ~ 7deg^ and 
an exposure time of 30s yielding a 4 minute response for the 
same error box. The Zwicky Transient Facility (ZTF[3) will 
be an extension of the Palomar Transient Factory (PTF; 
7 deg^ FoV) with an increased 47 deg^ FoV and an or¬ 
der of magnitude improvement in survey speed. gotcEB, 
a dedicated optical follow-up instrument with a proposed 
initial configuration of 18 deg^ instantaneous FoV will re- 
quire an exposure time of around 5 mins to reach 21 mag 
dsteeghs &i Gallowavll20i5l l: however, for this instrument a 
shallower inventory of the large error ellipses may be possi¬ 
ble. Additionally, BlackGElVo, expected on-line from 2016, 
aims to cover 40 deg^ in a s eeing-limited configuration with 
an exposure time of 5 mins jGhosh fc Nelemansll201^ . 

Gonsidering the FoVs given in table [T] it is immedi¬ 
ately apparent that for fast response observations, unless 
a wide FoV instrument is on-source, one may require a 
larg er network than LHV t o decrease the localisation er¬ 
ror dManzotti fc Diet3l2012l ~) and thus, reduce the number 
of tiles required to image the full GW error box. 

In addition to the facilities outlined in Table |4l a 
number of other facilities are in development for GW 
follow-up and could be operational during the later era of 
aLIGO/AdV. One such facility is RAMSES which will con¬ 
sists of a set of 16 X 40cm telescopes and will be equipped 
with rapid cameras capable of readout times from 100ms to 
several minutes; it will be able to reach magnitudes of 15 in 
Is and mag 20 within around 10 mins. Each telescope has 
a 2.5 deg^ FoV, resulting in a combined 100 deg^ FoV; this 
in combination with a slew time of around 5s would enable 
instantaneous coverage of a large proportion of a typical 
aLIGO/AdV error region. RAMSES has received funding 
and will be based in National Aures Observatory in Alge- 
ri 43 with a prototype built in France. 


4.3 Radio follow-up 

Table [S] provides specifics for some of the radio instruments 
that have signed MoUs with LIGO/AdV. A period of signifi¬ 
cant investment in facilities and infrastructure and advances 
in high-speed computing over the last few years will enable 
this band to make a telling contribution to time-domain as¬ 
tronomy. The instruments included in this table have the 
potential to be leading transient detectors by virtue of their 
large FoV, high sensitivity and significantly for low-latency 
follow-ups, the capability to respond from 10s of seconds 
up to within a minute. The triggering methodology has al¬ 
ready been implemented by th e Arcminute Microkelvin Im¬ 
ager fAMI: IStaIev et ;rni2013l'l who have demonstrated typ¬ 
ical response times for follow-up of Swift-'BAT triggers of 
order 5 mins. 

The ASK AFf3 Survey fo r Variables and Slow Transients 
(VAST; iMurphvlt aP |2013|'I is a Survey Science program 
that will operate in Western Australia. Phased array feed 
technology on each of 36 antennas 12-m diameter will allow 


http://www.ptf.caltech.edu/ztf 

http://goto-observatory.org/ 

https://astro.ru.nl/blackgeni/ 

http://www.craag.dz/ _ 

Australian SKA Pathfinder lASKAP: Ijohnston et ^l200it) 


an instantaneous field of view of 30 deg at 1.4 GHz, with 
a resolution of 10-30 arc sec at 300 MHz bandwidth, and a 
frequency range of (0.7 - 1.8) GHz. The VAST transient de¬ 
tection pipeline will operate on an imaging cadence of 5-10 
seconds at its fastest down to a cadence of minutes depend¬ 
ing on available supercomputing resources. Longer cadences 
of hours to months can also be conducted through repeated 
observations of selected fields. When commissioning is com¬ 
plete, ToO response timescales could be as short as minutes 
when operating in fully automated mode. This accompanied 
by the large FoV could provide good follow-up opportunities 
for very low-latency triggers. _ 

In the Netherlan ds, Apertif (e.g. IVerheiien et al.ll2008l : 
lOosterloo et al.lbOld ') will significantly increase the field of 
view of the Westerbork Synthesis Radio Telescope and will 
operate at 1.4 GHz. In standard observing modes, the field 
of view will be 8 deg^, but has the potential to be exten ded 
to 100 deg^ in a “fly’s eye” mode dOosterloo et al.l 2010h. 

The Murchison Widefield Array fMWA: iTingav et al.l 
I 2 OI 3 II is a Western Australia based low-frequency radio tele¬ 
scope operating between 80 - 300 MHz. It consists of 2048 
dual-polarization dipole antennas arranged as 128 tiles. The 
array has no moving parts, and pointing is achieved by elec¬ 
tronically steerin a^^l . This provides the ability to retroac¬ 
tively point, at a reduced sensitivity, to a location on the 
sky in response to a GW event trigger. In addition, imple¬ 
mentation of a VOEvent triggering pipeline will reduce th e 
response of this instrument to 10s of seconds (lKaplanll2015l l. 

The MWA’s large FoV - 600 deg^ at 150 MHz defined 
by the FWHM of the primary beam - can be extended to 
1600 deg^ by taking the entire central lobe of the primary 
beam. At lower frequencies, the FoV scales proportionally 
to the square of the wavelength, so at lOOHz one can obtain 
a FoV ~ 3600, almost 10% of the celestial sphere at a single 
pointing llKaplan et al.l[^15l L Given a sufficient event rate, 
the MWA could be pointed at a single patch of sky guar¬ 
anteeing that 10% of GW triggers would be captured with 
zero latency. This improved FoV can be augmented by the 
MWA Voltage Gapture System (VCS flVemblav et al.l[2015h 
which would allow data to be simultaneously collected in 
a 2s cadence imaging mode and in continuous voltage data 
collection mode. The VGS can be used as a deep disk-based 
ring buffer with a 1 hr duration. On receipt of a GW trigger, 
voltage data around the on-source window could be siphoned 

off for post-trigger analysis. _ 

_ The Low Frequency Array iLOFAR: lvan Haarlem et al.l 

l2013ll is an interferometric array of dipole antenna stations 
distributed throughout the Netherlands and several coun¬ 
tries in Europe. LOFAR uses two distinct antennae types: 
Low Band Antennae (LBA) operating between 10 - 90 MHz 
and High Band Antennae (HBA) operating between 110 - 
250 MHz. The antennae are organised in aperture array sta¬ 
tions, with initial beam forming conducted on the station 
level. These stations have no moving parts and electronic 
beam-forming techniques allow for rapid re-pointing of the 
telescope. Due to strategies required to manage the high 
data rates from LOFAR, the FoV is constrained during the 
initial processing stages; it is highly dependent upon the 


Telescope functions including pointing are performed by elec¬ 
tronic manipulation of dipole signals. 
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station configurations chosen and the observing frequency. 
Larger FoVs can be attained by using shorter baseline con¬ 
figurations but with a loss of sensitivity and resolution. The 
maximum FoV (whole observable sky when operating with 
the LBA) will be monitored in a 24-7 commensal observing 
mode by (AARTFAACj3 (the Amsterdam-ASTRON Ra- 
dio T ransients Facility and Analysis Centre; IPrasad et al.l 
l2014h At the time of writing, LOFAR still requires man¬ 
ual intervention for external triggers. However, the expecta¬ 
tion is that VOEvent triggering will be enabled in time for 
the first LIGO science runs, reduce the latency to less than 
a minute. In addition, LOFAR are also developing rapid 
response strategies based around transient buffer boards 
llvan Haarlem et al.ll201a ') which, similar to the VCS used 
by MWA, can store small amount of the full observed data 
for later analysis. 

In the longer term, the Square Kilometre Arra^F^ 
(SKA) will be provide large fields of view and increased sen¬ 
sitivity, with construction scheduled to start in 2018. SKA 
phase 1 will comprise of two instruments SKA-low, operat¬ 
ing at 50 - 350 MHz, and SKA-mid, operating at 350 MHz 
- 14 GHz. 


5 A REVIEW OF EM EMISSIONS FROM BNS 

MERGERS 

A number of associated EM counterparts have been pro¬ 
posed for BNS mergers, many of which have already been 
presented in Figure 2. We will begin this section by expand¬ 
ing on the multi-messenger pathways and end products pro¬ 
vided in Figure 2 which could produce observable EM coun¬ 
terparts for binary coalescences. We will discuss the various 
scientific breakthroughs that could be achieved in the ad¬ 
vanced detector era through low-latency follow-ups and dis¬ 
cuss the type of EM observational bands and facilities that 
could benefit the most from low-latency GW triggers. 

Although the scope is wide, the relative brightness of 
the associated EM signals is dependent on poorly under¬ 
stood merger physics, in particular that of the post-merger 
central remnant object. Therefore, the prospects of detect¬ 
ing many of these EM counterparts are quite speculative. To 
address some of this uncertainty, in the latter parts of this 
section we use the best available EM data to indicate the 
possible flux of the possible EM counterparts. We do this in 
the context of EM telescopes already discussed in section 4 
which will be available for follow-up observations during the 
aLIGO era. 

Short GRBs have long been propo sed to originate 
from mergers of compact object binaries ( Pacz yr^M 19 ot; 
lEichler et al.l[l989bl : IPaczvnskil[l99ll : lNaravan et ah 1992bll . 
Recent broad-band observations have shown mixed host 
galaxy types, non-detection of supernova associatio ns and 
offset s of GRB locations from their host galaxies fsee lBerg^ 
|2014 for a review). These findings lend support to compact 
binary mergers as the progenitors for SGRBs. 

Generally the merger of BNSs can produce a number of 
different products; these can be divided into four categories: 


AARTFAAC will have a FoV of 10000 deg^ covering 25% of 
the sky. See http://www.aartfaac.org/ 
https://www.skatelescope.org/ 


1) A black hole. 

2) A hyper-massive NS supported by differential rotation; 
this is expected to collapse to a BH within lO-lOOrns. 

3) A stable NS that does not collapse. 

4) A supra-massive NS supported by rigid rotation; this 
is expected to collapse to a BH within 10-10000s with a 
typical timescale of hundreds of seconds. 


If the merger product is of the first two catagories: 
1) a black hole or 2) a temporal hyper-massive neutron 
star, the EM signals could include the SGRB and its af- 


terglow emission (lEichler et alJll989b 

R0SSW02: et al.ll201S: 

iGehrels et alJ l2005l: Barthelmv et al 

200 , 5l: Berger 1201 11 

and a long-lasting radio afterglow of the ejecta as it in- 

teracts with the ambient medium dNakar & PiranI I 2 OIII: 

iMetzeer & Bereeill2012l: IPiran et alj|2 
During a NS-NS merger or NS-BL 
tion of materials are tidally ripped anc 
These ejecta are mostly composed of 
which would heat the ejecta via deca 
powering an UV/optical/IR transient ( 

merger, a small frac- 
dynamically ejected, 
radioactive elements, 
V of r-process nuclei. 

Li & Paczvnskilll99^ 


iKulkarrii l2005l : iMetzger et al.ll20ld : iBarnes fc KasenI 12013 ) 


that has been termed a kilonova. Recent optical and near- 
infrared observations of a kilonova associated with GRB 
130603B have provided support for this scenario and thus a 
compact object m erger origin of SGRBs llTanvir et al.l[201^ : 
iBerger et al.ll201^ 'l. 

It is important to note that the SGRB and its after¬ 
glow are relativistic an d collimated dBurrows et al. I l2006l : 
IPe Pasauale et H] l2010l L so only detectable in preferred 
directions. Alternatively, the kilonova ejecta and its radio 
afterglow are essentially isotropic and also non-relativistic 
(due to the heavy mass loading) and therefore can be de¬ 
tected from any direction if the flux is high enough. 

The latter two categories, 3) and 4) can both lead to a 
magnetar as the merger product; this could occur through 
a small enough total NS mass in combination with a stiff 
equation of state of nuclear matter. This could be a stable 
magnetar (for case 3) or a supra-massive millisecond mag- 
neter (case 4), which would ver y likely be strongly magne- 


20061: iFan & Xul 20061: iGao & FanI 

20061: 1 Yu & Huang 2007: 

Zhang 2 OI 3 I: Giacomazzo & Perna 

20131: Hotokezaka et al. 

20 ia 

: Shibata et al. 20141: Laskv et al.l 20141: Ravi & Laskv 


2014fl . 


Plateaus in the X-ray afterglow light curves of about 
60% short GRBs have been attributed to such m i llisec¬ 
ond magnetar central engines dRowIinson et al]l2010al . l2013l : 
ILu et aI.ll2015l L The plateau emission of these SGRBs is in¬ 
terpreted as internal dissipation of the magnetar wind due 
to the dipole spindown of the magnetar (IZhang Sz Meszarod 
I2OOII : iTroia et al.ll2007l : IMetzger et al.ll 2011 ^ . The magnetic 
spin-down of the NS remnant provides an additional source 
of sustained energy injection, which greatly enriches the 
EM signals. A jet may be still launch ed via accretion from 
a disk to power th e SGRB emissio n jMetzgeiLgLalJ 2004 
Dessart et al. I I2009I : iLee et ahl l2009l : iFernandez fc Metzgerl 


2013tl . In the SGRB afterglow phase, some observed features 


including the extended prompt emission ((Norris_&_^onn^ 


2001 

: IMetzger et al.l 

20081) 

2005 

:lGamDana et al. 

to 

0 

0 

QL 


X-ray flares ( Barthelmv et al.l 
and more importantly, plateaus 
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observed in the X-ray are consistent with being powered by 
a millisecond magnetar central engine. 

Outside the solid angle extended by the SGRB jets, sev¬ 
eral potentially observable astrophysical processes have been 
predicted. First, self-dissipation of the magnet ar wind coul d 
power a bright early isotropic X-ray afterglow dZhan iEH). 
Next, heating by the magnetar wind (expected to be more 
significant than the radioactive decay) co uld significantly 
enhance the kilonova emission component (IYu et al.ll2013l : 
iMetzger fc Pirdl20l4l . Finally, the magnetar wind could ac¬ 
celerate the ejecta to a mildly or even moderately relativistic 
speed, giving rise to a strong broad-band after glow emis¬ 
sion upon interaction with the ambient medium dOao et al.l 
l2013h . 

Finally, a remarkable new class of extremely bright and 
short-duration (millisecond timescale) transients called Fast 
Radio Bursts (FRBs) could be a product of BNS merger 
dThornton et al.ll2013l : [borimer et al.ll2013h . An FRB could 
be emitted through a BNS merger in two possible ways: 

(1) during the BNS merger a coherent radio emission could 
be produced through magnetic braking as the magnetic 
fields of the NSs are synchronize d in binary rotation at the 
time of coalescence l|Totanill2013l l: 

(2) if a supramassive NS is the merger product, it would 
collapse into a black hole after loosing the centrifugal sup¬ 
port and emit an FRB dZhand l2014h due to the ejection 
of the magnetosphere a s it implodes into the black hole 
(|Falcke fc Rezzollall201-j l. 

The latter scenario is evidenced by an abrupt drop in 
the X-ray plateau emissions of some long and short GRBs. 
Within this latter scenario, an FRB along the jet direction 
is expected after an BNS merger event that forms the supra¬ 
massive NS. 

In the remainder of this section, we will discuss each 
of the potential electromagnetic counterparts in more depth 
and compare the expected emission to the capabilities of the 
instruments described in Section 4. First, in Section 5.1, we 
outline the methods used to quantify the typical emission 
we might expect. Section 5.2 considers the on-axis emission 
expected from our understanding of SGRBs and Section 5.3 
discusses the isotropic kilonova counterpart. In Section 5.4 
we describe the impact of changing the central engine formed 
by the BNS merger from a black hole to a magnetar on 
the standard emission as discussed in Sections 5.2 and 5.3. 
Finally, in Section 5.5, we consider the possibility of early- 
time FRBs associated with a BNS merger and how low- 
latency triggers and existing instrumentation can constrain 
these predictions. 

6 ASSESSING PROMPT EM FOLLOW-UPS 

FROM LOW-LATENCY GW TRIGGERS 

6.1 Extrapolation of fluxes for observed EM 
counterparts 

In this section we explore the benefits of low-latency GW 
follow-up observations by extrapolating observed emissions 
from known sources associated with BNS mergers to the 
observational range of ground based GW detectors. SGRBs 
are the most likely observed population to be associated 
with BNS. Therefore, we can estimate the typical fluxes that 


could be expected at different energies by extrapolating from 
known SGRBs associated with BNS at a distance of 200 
Mpc. The results are shown in Table [S] however, due to the 
nature of the SGRB observations, these numbers focus on 
the on-axis emissions. These values are representative of the 
expected on-axis signals but we do not give the full range 
of expected fluxes; this requires detailed modelling, beyond 
the scope of this paper. Where possible, we utilise GRB 
130603B as it is the only unambiguously SGRB associated 
with a BNS merg er with a spectroscopic redshift from the 
optical afterglow jBarthelmv et al.l l2013l : iFong et al.l l2014l : 
Ide Ugarte Postigo et al.ll2014^ . GRB 130603B does not ex¬ 
hibit the all of the emission combinations we consider and, 
in those cases, we use other SGRBs with reasonable host 
galaxy redshift constraints. 

The X-ray, 7-ray and radio flux estimates are obtained 
by extrapolating to the required observation time from bro¬ 
ken power law fits to the lightcurves of the GRBs, which are 
converted into rest-frame luminos ities using a k-corr ection 
and the observed spectral slope (iBloom et al.l I2001I ') . The 
observed flux at a distance of 200 Mpc is calculated by 
converting the rest-frame luminosities using an inverted k- 
correction for the new distance. The stan dard 7-ray emis¬ 
sion is extrapolated from GRB 130603B dBarthelmv et al] 
l2013tl . while the long lasting VHE 7-ray emission is ex- 
trapolated from Fermi L AT observations of GRB 090510 
dDe Pasauale et ^ l2010l l. In Sections 5.2.2 and 5.4, we 
consider three different X-ray counterparts and use three 
SGRBs with reasonable redshift constraints to quantify the 
typical fluxes given in Table [G] namely: GRB 050509B for 
a SGRB without energy injection in the X-ray afterglow 
jGehrels et akll^OOSl : iTunniciiffe et al.l[20l3l . GRB 130603B 
for a long lived afterglow showing a clear signature of energy 
injection associa t ed with a stable magnetar central engine 
dFan et al.l l2013l : Ide Ugarte Postigo et al.l I20l4 iFong et ^ 
l2014h and GRB 080905A for an X-ray lightcurve with 
the signature of energy injection from a short-lived cen- 
tral engine, such as an unsta ble magnetar central engine 
dRowlinson et all l2010bl . l2013f l. The radio fluxes ^iven in 
Table [S] are extrapolated from GRB 130603B llFong et al.l 

I2OI4II . 

The optical magnitudes are determined at 200 Mpc us¬ 
ing the standard formalisrrF^. The optical afterglows are 
considered in Section 5.2.2; we typically observe very faint 
and rapidly fading optical afterglows to SGRBs and we rep¬ 
resent these counterparts in Table Blby extrapolating the op - 
tical counterpart to GRB 080905A ^Rowlinson et al]2010bl) . 
GRB 130603B had an unusually bright optical counterpart 
and we use this GRB to represent the bright optical after¬ 
glow in Table (6] Additionally, in Section 5.3, we consider 
the isotropic kilonova counterpart to GRB 130603B which 
would have an observed magnitude of ~21 in the J band, 
but undetected in the R-band to >23.7, at a distance of 200 
Mpc. 

In the following Sections, by combining the results from 
Tables OEl we discuss some possible observational opportu- 


The standard formula is, m — M = 51ogT> — 5, with m the 
apparent magnitude, M is the absolute magnitude and D is the 
luminosity distance in parsecs 
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Table 6 . The typical expected emission from an on-axis compact binary merger at a distance of 200 Mpc for given observing latencies 
(- = no emission expected, ? = possible emission, refer to text for more details). We provide different predictions for th e X-ray emission 
repre senting the range of different lightcurves observed from SGRBs (interpreted as different energy injection signatures: iRowIinson et al.l 
l2013h and separate the SGRBs with bright and faint optical counterparts. 


Wavelength 


EM follow-up time 



0 s 

2 min 

10 min 

1 hour 

9 hours 

VHE 7 -ray (>100 MeV, ph cm“^s“^) 

1.1 X 10^ 

1.2 X 10“^ 

< 1.7 X 10“® 

- 

- 

7 -ray (15 — 150 keV, ph cm~^s“^) 

3.2 X 10^ 

7 

7 

- 

- 

X-ray - no energy injection (0.3 - 10 keV, erg cm~^s“^) 

- 

7 X 10-12 

4 X 10-12 

2 X 10-12 

9 X 10-1® 

X-ray - unstable magnetar (0.3 - 10 keV, erg cm~^s“^) 

- 

3 X lO-'l 

6 X 10-11 

5 X 10-1® 

1 X 10-1® 

X-ray - stable magnetar (0.3 - 10 keV, erg cm~^s~^) 

- 

1 X 10“® 

6 X 10-2 

2 X 10-2 

8 X 10-11 

Optical - bright (R band, mag) 

7 

17.0 

16.8 

16.5 

16.9 

Optical - faint (R band, mag) 

7 

17.9 

19.0 

20.2 

21.8 

Radio (6.7 GHz, mjy) 

7 

7 

7 

7 
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nities and implications provided through a low-latency win¬ 
dow. 


6.2 Short duration GRBs 

6.2.1 Probing the progenitor models of Short GRBs 


Coincident EM/GW observations of SGRBs can prove with¬ 
out doubt that BNS mergers are the progenitors behind 
these bursts. It is quite possible that 1-2 coincident events a 
year may be acc essible at full aLIGO/AdV sensitivity using 
GRB satellites ([Metzpr fc Bergei]l2012l : ISiellez et ai.ll2014f 
iRegimbau et ahl 20151) . Such a verification could be pro¬ 
vided by a wide FoV instrument such as Fermi. Low-latency 
follow-ups could also enable multi-wavelength measurements 
of the prompt emission allowing the underlyin g central en¬ 
gine and emission mechanism to be probed llElliott et al.l 

I2OI4I) . 

A particular obstacle faced by Swift-led observations of 
SGRBs is that the observed 7 -ray emissions typically have 
observed durations shorter than 2 seconds. Therefore even a 
GW trigger with a latency as short as 10 s wold not be suf¬ 
ficient to allow optical telescopes to recover emission during 
the prompt phase. However, a GW trigger sent out during 
the inspiral phase would allow multi-wavelength telescopes 
be on target before the SGRB. Such a scenario could poten¬ 
tially provide a breakthrough in our understanding SGRB 
prompt emission physics. 

In fact a number of foreseeable opportunities to study 
SGRB physics exist that could be plundered through low- 
latency follow-ups. Within the GRB theoretical framework, 
panchromatic emission during the prompt phase and early 
afterglow phase is predicted, including the emission from 
the prompt emission site through internal shocks (e.g. 
iMeszaros fc Reegll993h . Indeed, for long duration GRBs, 
panchromatic e mission has been observed in a ha ndful of 
Swift GRBs (e.g. lPaee et al.l[2007l : IZheng et al.ll2012l) thanks 
to their long durations or the existence of a precursor 
that triggered the Swift BAT. In particular, prompt opti¬ 
cal flashes have been observed in a handful of long GRBs: 
e.g. the Naked eye burst GRB 080319B dRacusin et al.l 
I 2 OO 8 I). GRB 990123 llAkerlofet all I 1999 I) . GRB 041219A 
JVestrand et al.ll2005l) . and GRB 130427B (IVestrand et al.l 
I 2 OI 4 I ). However. so far no emission outside the 7 -ray band 


has been detected during the prompt emission phase of 
SGRBs. Low-latency GW triggers 10s of seconds before 
merger could make this possible for the first time. As noted 
in section 4.2, such observations will become more feasible 
once the error region is decreased through an expanded GW 
network in the latter part of the advanced GW detector era. 
For a particularly bright or close counterpart, the wide FoV 
and fast response of an instrument such as GWAC could 
prove to be highly valuable. 


The high energy reach of Fermi has allowed access to 
photons with energies in access of 10 GeV in both long and 
SGRBs. The fact that Fermi LAT discovered a photon of 
energy 31 GeV during the prompt phase of SGRB 090510 
dAckermann et al]l20ld ) is quite significant for low-latency 
follow-ups at VHE. In TablelS] we provide the expected VHE 
(> 100 GeV) 7 -ray flux and 7 -ray flux for a putative SGRB 
at 200 Mpc coincident with a GW source. A low latency 
trigger GW trigger would alert VHE and 7 -ray detectors to 
search for prompt EM counterparts. In the 7 -ray band the 
wide FoV of an instrument such as Fermi would be advanta¬ 
geous for such rapid follow-ups. At higher energies, instru¬ 
ments such as CTA coul d provide follow-ups in wide-held 
mode dBartos et al.ll20l3 ): such observations would be very 
rich in high energy photons and could provide v aluable clues 
on the prompt emission mechanisms of SGRBs dinoue et al.l 
I2OI3I : iBartos et al.ll20l3 ) . 


Additionally, we note that there are a small popula¬ 
tion of SGRBs that exhibit extended 7 -ray emission for 
hundreds of seconds fo llowing the initial short burst (e.g. 
iNorris &: Bonneiill2006l ). This sub-population of SGRBs are 
also believed to originate from compact binary mergers, and 
hence we might expect to observe this extended emission 
in some of the triggers even if the bright short-hard emis¬ 
sion is missed by the 7 -ray detectors due to the inadequate 
time for slew from the GW trigger. The possible response 
times provided for CTA in Table 3 suggest that VHE follow 
ups of these emissions are quite possible. Such observations 
could po tentially unravel the exact nature of this e xtended 
activity dinoue et al.ll2013l : IVeres fc Meszar^l2014l ). 
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6.2.2 Capturing the afterglow and late-time central engine 
activity of SGRBs 

By comparing the expected X-ray counterparts given in Ta¬ 
ble [6] to the instruments described in table |3] we show that 
existing X-ray instruments such as Swift XRT or SVOM 
MXT could detect the early X-ray afterglow for the SGRBs 
even without energy injection. Unfortunately, if the X-ray 
observations start with a latency longer than ~10 minutes 
following the BNS merger, or they require a long time to 
tile out regions of the sky due to having a small FoV, we 
are likely to miss the X-ray counterparts for SGRBs with¬ 
out prolonged energy injection in their lightcurves (such as 
those described in Section 5.4). 

We consider both bright and faint optical afterglows in 
Table [6] and note all of the optical instruments would be 
capable of observing the earliest onset (or the rise) of the 
optical afterglow assuming the 17th magnitude extrapola¬ 
tion based on a burst like GRB 130603B. However, GRB 
130603B was unusually bright at early times, with a likely 
optical plateau phase, so the optical emission may have been 
contaminated with ener^ inje ction from the central engine 
dde Ugarte Postigo et al. 120141 '). Hence, these values are op¬ 
timistic limits. Most SGRBs typically have a fainter counter¬ 
part, such as GRB 080905A, so will require deeper optical 
observations. The short optical reach and short exposure 
times of Pan-STARRS and ZTF suggest these instruments 
could be valuable in such endeavors. The larger FoV of ZTF 
would be more suited to searching the error regions from ear¬ 
lier GW triggers and could facilitate a co-ordinated search 
strategy with an instrument such as Pan-STARRS capable 
of deeper searches but more limited by FoV. 

In the radio domain, there have been 3 detections of 
radio emission associated with SGRBs (GRB 050724,GRB 
051221A and GRB 130603B) as well as a forth case, GRB 
120804A, due to host galaxy emission. In Table[6l we provide 
fluxes based on the radio counterpart of GRB 130603B. All 
of these detections have been at late times (> 7 hours) in the 
high frequency domain (> 4.9 GHz) with fluxes around the 
mjy level. These fluxes could be observed by instruments 
such as ASKAP or APERTIF. The peak in the radio after¬ 
glow would be expected to occur on typical timescales of 
days - years. At low radio frequencies, we may also expect 
to observe the afterglow on months - years timescales using 
facilities such as MWA and LOFAR although this emission 
has not been observed to date for a SGRB. 

During the early afterglow phase, a short-lived back¬ 
ward propagating reverse shock component can occur once 
the shocked matter encounters the interstellar medium. This 
may dominate the optical and radio band under certain con¬ 
ditions, e.g. if the reverse shoc k region is more m agnetized 
than the forward shock region (IZhang et al.ll2003f ). 

On short timescales, an early radio flash (< 1 day) 
originating from a reverse shock may be expected. Reverse 
shocks have been observed for a number of long GRBs, 
peaking a few hours after the G RB, (e .g. GRB 130427A 
lAnderson et al.l l2014l : IPerlev et al.l l2014ll and may already 
have been observed in some SGRB radio afterglows (e.g. 
a radio excess was potentially observe d at early times for 
GRB 130603B; iGonipertz et al.ll2015a^ . Early radio obser¬ 
vations, within the first hours of the BNS merger, cap¬ 
turing the rise of the reverse shock at different radio fre¬ 


quencies would allow detailed modeling of this mechanism 
and would be a valuable probe of the jet models, macro¬ 
physics within the jet and the surrounding environment of 
the GRB (as conduc ted for the long GRB 130427A by, e.g., 
Ivan der Horst et al. lliH). With its expected rapid trigger¬ 
ing capabilities, ASKAP has the capability to detect an early 
reverse shock. 


6.3 Kilonova 


The peak frequency and brightness of a kilonova criti¬ 
cally depend on the amount of ejected materials, and more 
importantly, the opacity of the eject a. Early calculations 
adopted a relatively s mall opacity (e.g . Li fc Paczvhskiiri998l : 
iMetzger et al.l I2OI0I : iBauswein et al.[ l2013l L and predicted 
an optical transient with a peak flux about 1000 times of 
nova luminosity (and hence, the term “kilonov a” is coined) 
d Metzger et al. I2OIC1I L iBarnes &: KaserJ (I2OI3I I pointed out 
that the existence of heavy elements, especially lanthanides, 
greatly increase the opacity, so that the kilonova signature 
would peak at a lat er time and in a softer frequ ency, typi¬ 
cally in IR (see also lTanaka fc Hotokezal^l2013l L 

The properties of the kilonova depend on the proper- 
ties of the post-merger compact object. The calculations of 
iBarnes fc KasenI (I 2 OI 3 I I are relevant for the case of prompt 
formation of a black hole. If the merger product is a supra- 
massive neutron star (sustained by rigid rotation for 100s 
of seconds or even longer, SMNS) or even a hyper-massive 
neutron star (sustained by differential rotation for 100s of 
milliseconds, HMNS), neutrinos from the NS would raise the 
electron fraction and reduce t he opacity, so that the kilonova 
signature is bluer and earlier (iMetzger et al.ll2015l L 

As discussed earlier, an observatio nal candidate for a 
kilonova was claimed for SGRB 130603B llTanvir et al.ll20r^ : 

I Berger et aklfioi^ l with one HST excess data point in the 
near IR band 9 days after the SGRB. As shown in Section 
5.1, this kilonova is equivalent to a J band magnitude of 21 
at 200 Mpc and an R band magnitude of >23.7. None of 
the optical telescopes given in Table 4 have near IR capa¬ 
bilities and are not expected to detect an R-band counter¬ 
part at this observation time. Detailed modelling will be re¬ 
quired to estimate the R band fluxes and peak timescales for 
these instruments (see the w ork by iTanaka fc Hotokezakal 
I2OI3I : IBarnes &: KaseiJ l2013l L Besides this one case, two 
other candidates have been reported. One is an excess 
around 13.6 days iri the HST F814W-band following GRB 
060614 dYang et al.ll201a '): another is a significant optical 
rebrightening feature around 1 day following GRB 080503 
(iPerlev et al]l200^ : iGao et al.ll201^ ~). More data are needed 
to reveal the diverse kilonova phenomenology. A low latency 
trigger would facilitate broad-band follow up observations to 
catch the kilonova signals, especially the early ones due to 
the small opacity in the ejecta. This will be a challenge due 
to the lack of sensitive instruments in the NI R band with 
FoVs large enough to survey the error regions (iBartos et al.l 
I2OI5II . As discussed in the next section, an associated X- 
ray emission could be possible; this would allow a wide field 
instrument such as ISS-Lobster to narrow down the search 
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6.4 Magnetar formation 

The X-ray plateaus attributed to a millisecond magne¬ 
tar central engine are observed up to 100-10000s after the 
bursts. The estimates in Table[^show that if a stable magne¬ 
tar is produced, GW triggers with a moderate latency could 
easily allow an X-ray instrument such as Swift to detect 
the plateau emission. However, should an unstable magne¬ 
tar be produced, a GW trigger with a 40s latency would be 
required enable an EM follow up within a few 100s. Such 
observations of the platau phase in combination with GW 
information c ould provide important clues of t he NS equa¬ 
tion of state (|Laskv et al.ll2014 : iLii et al.ll201!Th . 

In addition to the X-ray plateaus. X-ray flares observed 
in the first few hundred seconds following short GRBs have 
also been interpreted as resulting from magnetic activity 
in a post merg er differentially rotating supramassive NS 
dPai et al.ll2006l l. Prompt follow ups of the GW signals from 
NS mergers would provide significant evidence for this for¬ 
mation channel of magnetars. 

In the massive proto-magnetar central engine model, 
the potentially wide-beamed magnetar win d is expected 
to produce a bright broadband afterglow (IZhand |2013l : 
iGao et ahllioT^ . This could result in an early afterglow in 
the optical (possibly as bright as 17th magnitude in R band) 
and X-ray bands (up to 10~^ — 1 0~^ erg s“^ cm“^) for a 
source at ~ 300 Mpc dZhan 3 12013) . 

The existence of a supra-massive NS can also supply 
continuous heating to the ejecta which can, depending on 
the spindown parameters of the neutron star, greatly en¬ 
hanc e the kilonova signal (|Yu et al.l I 2 OI 3 I : iMetzger fc Pirol 
l2014h . The detectability also depends on the observational 
geometry. For an on-axis binary merger, the signal may be 
fainter than the SGRB optical afterglow. However, for off- 
axis mergers, this would be the main UV/optical/IR sig¬ 
nal from the event. The radio emission is also expected to 
peak earlier and brighter in this central e n gine model (e.g. 
Gao et al. I I2OI3I : IMetzger &: Bower! I20l4 iGompertz et al.l 
2015bl ). TTese predicted optical and radio counterparts are 
potentially within the sensitivity and field of view capabili¬ 
ties of many of the optical and radio telescopes described in 
Tables 4 and 5. 

Additionally, recent suggestions that the observed 
prompt emission does not necessarily signal the start of 
the X-ray emission could ha ve important cons eq uence s 
for low latency fo llow-ups. iRezzolIa fc Kurnarl (l2014l ): 
ICiolfi fc Siegell l|2014h have postulated that the X-ray emis¬ 
sion could start before the observed short hard spike emis¬ 
sion of short GRBs. These models predict that X-ray emis¬ 
sion from a proto-magnetar could occur 100-1000s prior to 
the prompt GRB. 

It is clear that low latency observations will play a ma¬ 
jor part in testing these models, but these will require more 
rapid triggering pipelines than those achievable by facilities 
such as Swift and SVOM-MXT (see Table 3). The employ¬ 
ment of a wide FoV X-ray instru ment with rapid p ointing 
capabilities such as ISS-Lobster (jCamp et al.ll2013l ) would 
be very timely for the advanced GW detector era. 


6.5 Fast Radio Bursts 

The possibility of associations with FRBs and SGRBs 
and/or BNS merger s is still highly speculative F3 . 
iThornton et al.l (l2013l ) has suggested that the event rate 
of FRBs (~ 10® gal“®yr“®) is inconsistent with GRBs and 
compact star mergers, but could be in agreement w ith soft 
gamm a-ray repeaters or core-collapse supernovae. iTotanil 
1I2OI3I) argued that including the effects of cosmological 
time dilation and merger rate evolution would make the 
FRB event rate in agr eement with the o ptimistic double 
NS merger event rate (lAbadie et ahllioiol ). One caveat is 
that these models suggest an essentially isotropic emission of 
FRBs. However, due to the existence of the r-process ejecta 
launched before the coalescence, the FRB signal may not be 
able to escape the m erger environment in most of the solid 
angles. IZhand (|2014| ) acknowledged this and propose that 
the GRB- or merger-associated FRBs only account for a 
small fraction of the entire FRB population, and only those 
mergers whose jet direction coincides the observer’s viewing 
direction (i.e. the short GRB direction) would be associated 
with FRBs. The majority of FRBs are attributed to those 
SMNSs whose delay time scale to collapse into a black hole is 
much lon ger, and therefore not ass ociated with GRBs or su¬ 
pernovae (|Falcke fc Rezzollall2014h . This is certainly an area 
where investments in GW low-latency trigger and prompt 
radio follow-up observations could provide significant scien¬ 
tific insight. 

As discussed earlier, a particular obstacle for short du¬ 
ration events of order a few seconds, is that even a latency 
as short as 10 s would not be sufficient to alert a 7 -ray in¬ 
strument in time to recover the prompt emission. This how¬ 
ever can be circumvented in the radio band, by electroni¬ 
cally steered radio arrays such as LOFAR, MWA and, in 
the future, SKA-low. An instrument such as MWA could 
cover the entire GW error region in less than a minute of re¬ 
ceiving a trigger. Although LOFAR typically has a smaller 
field of view than MWA (this is dependent upon observ¬ 
ing configuration and frequency), it is also capable of tiling 
out large regions of the GW error box using multiple si¬ 
multaneous pointing directions and can conduct commensal 
beam-formed observations to probe higher time resolutions. 
In addition, as MWA and LOFAR operate at low frequen¬ 
cies the signal will be delayed by propagation through the 
ionised intergalactic/interstellar media. Assuming a disper¬ 
sion measure of 1000 pc cm®, the delay could be of order 3 
minutes at 150 MHz - sufficient time to point the telescope at 
the error region. With the voltage capture system on MWA 
and the transient buffer boards on LOFAR, in combination 
with the low latency triggers and signal delay time, it will 
be possible for data to be obtained that can constrain emis¬ 
sion from seconds to minutes prior to the merger. SKA-low 
is expected to have many of the capabilities of LOFAR and 
MWA but with a much higher sensitivity, making it the ideal 
future instrument for these triggers. For higher frequency in¬ 
struments such as ASKAP and APERTIF, the delay would 
be significantly shorter - of order 8 s at 700 MHz - requiring 
much more rapid response times. Once on target, these ra- 

However, some models for coherent radio emission have been 
proposed the late stages of a compact binary i nspiral phase or 
during their coalescence (e.g. lUsov fc Katal2000ll . 
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dio telescopes would ideally continuously monitor the entire 
field for several hours as the different predicted FRB mech¬ 
anisms cover a range of timescales, both a few seconds prior 
to the merger and hours following it, and these observations 
would also be able to constrain the occurance of other types 
of coherent radio emission. 


7 DISCUSSION AND CONCLUSIONS 

The first gravitational wave observation has provided a path¬ 
way to the highest energy density events in the Universe and 
transform our understanding of the behaviour of matter and 
space-time under the most extreme conditions. Coalescing 
systems of NSs are widely expected to be the progenitors of 
short GRBs. Early multiwavelength follow-up observations 
of these sources would allow one to probe and unravel the 
mechanisms behind these events beyond the GW regime; the 
potential for breakthough science is without doubt. 

In this study we have simulated the expected localisa¬ 
tion error regions corresponding to different times before the 
merger of a coalescing system of NSs. To examine the most 
probable low-latency EM observations we have considered 
the benefit provided by: a) larger networks of GW detec¬ 
tors, corresponding with improved localisation; b) improved 
GW triggering latencies for sending out GW triggers. 

The Tables shown in Figures [5] and [6] illustrate some 
of the multi-messenger observations that could be achieved 
through low latency GW triggers assuming a detection 40s 
before merger. We note here that these tables may well be 
riddled with caveats due to the different operating scenarios 
and evolution of the different wavelength instruments. Fur¬ 
thermore, it would require very detailed modelling to predict 
exact fluxes of the different EM emission mechanisms; well 
beyond the scope of this work. The goal of this study is 
to highlight the type of breakthrough science that could be 
realised through the low-latency domain and, using a wide 
selection of EM instruments with agreements in place with 
aLIGO/AdV for follow-up observations, determine the most 
likely avenues for success. 

FigureOpresents a conservative scenario; that of a LHV 
network corresponding with a 1000 deg^ error region (as¬ 
suming the 50% case provided in Table [2]) and a 40s GW 
triggering latency to send out the alert. At low frequency 
radio and both high and low-energy y-ray, the large FoV 
and fast response potential show that post-prompt (> 5s) 
breakthrough observations could be realised. 

In the optical, even for Pan-STARRS and ZTF, both 
large FoV instrument with fast responses, to cover the 1000 
deg^ error region with an exposure times of 30s would re¬ 
quire > (4000;600)s respectively. For longer latency follow¬ 
ups, optical instruments available in the latter stages of 
the aLIGO/AdV era include BlackGEM which has a 40 
deg^ FoV and a depth of 22 mag and a possible fu¬ 
ture configuration of GOTO, which could employ a 40- 
50 deg^ instantaneous FoV and sensitivity to 20-21.5 mag 
dSteeghs &: Gallowavll201^ '). 

To catch the early prompt stage, assuming a a 40s GW 
triggering latency, even instruments with a large FoV will 
be challenged; this scenario certainly provides motivation for 
speeding up the trigger pipelines of both the GW observato¬ 
ries as well as the high energy satellites. We note that for the 


low-frequency radio, one bonus in regards to EM follow-ups 
is that the dispersion of coherent radio signals can introduce 
a significant delay of up to minutes (as discussed in Section 
16.511 : as the existence of early time counterparts to coalesing 
BNSs in the radio band are still highly speculative, we have 
taken a conservative approach and neglected this effect, as¬ 
suming only the physical time in our estimates. 

A more optimistic scenario is presented in Figure [6] 
which shows how breakthrough science can be achieved with 
an almost instantaneous Is GW triggering latency and a 
fully expanded advanced detector network, LHVJIA. For 
the 182 deg^ error region the post-prompt (> 5s) situation 
is much more optimistic in the optical and high-frequency 
radio bands due to the FoVs of radio instruments with re¬ 
sponses <mins and wide FoV telescopes such as ZTF that 
could survey this error region in around 2 mins. In the X- 
ray, small FoV instruments such as Swift are hindered by 
the large GW error region. Although considerable efforts 
have been undertaken to remedy this situation through the 
use a requirement for rec onstructed sky loc ations to over¬ 
lap with nearby galaxies dEvans et al.ll20f3 l. this coupled 
with the slow response of high-energy satellites presents a 
significant bottleneck in a band for which many of the post 
merger EM emissions are known to exist. It emphasises the 
importance to have a wide EoV X-ray instrument capable of 
rapid pointing; one potential instrument that could remedy 
the situation is the NAS A proposed 30 deg ^ EoV soft X-ray 
instrument ISS-Lobster dGamp et al.ll20ll^ . 

Eor prompt (< 5s) observations the wide EoV of Eermi 
means this instrument could provide a coincident observa¬ 
tion with a GW trigger - this scenario is considered in the 
low gamma-ray ( < 10 MeV) region of Figures [5] and [6] How¬ 
ever, we note the slow response of other satellites operating 
at the I 0 W -7 will hinder prompt follow-ups. If such a low GW 
triggering latency could be achieved and instruments such 
as Swift/Fermi could employ fast triggering pipelines, and a 
wide FoV X-ray instrument could be employed, low-latency 
EM followups with a full GW global network would enjoy 
almost complete coverage for the models we have considered. 

We note that the scenarios we have outlined in Figures 
[S] andassume the GW error regions 40s before merger; this 
choice was to provide the EM facilities with the maximum 
time for response. However, as shown in Tabled closer to the 
merger both the error regions contract and the percentage 
of detections sources increases. Eor sources detected early in 
the inspiraling stage, subsequent alerts would most likely be 
sent out as the error box decreases in size; this suggests that 
a co-ordinated observational and alerting strategy involving 
both the GW network and both wide and narrower FoV in¬ 
struments would be optimal procedure for maximising GW 
detections. 

Our results support the argument that much break¬ 
through science can be achieved though a fast response. This 
has been evidenced in the post-Swift era with prompt follow¬ 
up of GRB triggers producing such unexpected phenomena 
as flares, plateaus, optical flashes and reverse shocks. His¬ 
tory therefore suggests that the lasting legacy from the era of 
prompt multi-wavelength GW follow-ups will undoubtedly 
be provided by the unexpected surprises we encounter. 
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SIGNAL 

TARGET 

TIME RANGE 

TYPE 


(s) 


SGRB Prompt 

0-5 

Jet Related 

SGRB Afterglow 

10 -> 


Reverse Shock 

60-10000 

Isotropic Ejecta 

Kilonova 

1000 -> 

Related 

Ejecta 

Afterglow 

10000 -> 


Magnetar Wind 
Dissipation 

1-1000 

Other Signals 

FRB (early) 

-5 -10 


FRB (late) 

10-10000 


Gamma 


High 

(> 10 MeV) 




■ 

Response fast enough 
FoV within range 


Fast response 

FoV too small 

✓ 

Observed 


Response too slow 

FoV within range 


Response too slow 
+ FoV too small 

7 

Predicted 


Figure 5. Some of the multi-messenger observations that could be achieved through low latency GW triggers. All scenarios assume a 
detection 40s before merger by an LHV network - corresponding with a 1000 deg^ error region - and a 40s GW triggering latency to send 
out the GW trigger. The signals are grouped by into jetted, isotropic or unknown; the approximate emission time ranges around merger 
are also given. As shown by the legend, the ticks indicate the EM observations with observational support; the question marks indicate 
those that have been proposed through theoretical models. All possible observations in the table are given a color code; good news is 
given by dark shaded cells that indicate there are EM instruments available in the band with sufficiently fast responses, exposure times 
and FoVs to cover or tile the GW error box within the emission time-range; at the other extreme, light shaded cells indicate that EM 
follow-up observations will be quite challenging. Between these two extremes there is scope for improvement; for example ToOs could 
potentially be speeded up or wider FoV instruments could come on line. 



Figure 6. As for Figure [5] but assuming an expanded LHVJIA network - corresponding with a 182 deg^ error region, and a Is GW 
triggering latency to send out the GW trigger. This figure shows what could be achieved by EM instruments given a ~40s head start to 
get on-source. 
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